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We evaluate the form factors governing the exclusive decays B —> pin, B —> a\(,v, by using an 
effective quark-meson lagrangian. The model is based on meson-quark interactions, and the compu¬ 
tation of the mesonic transition amplitudes is performed by considering diagrams with heavy mesons 
attached to loops containing heavy and light constituent quarks. This approach was successfully 
employed to compute the Isgur-Wise form factors and other hadronic observables for negative and 
positive parity heavy mesons and is presently used for exclusive heavy-to-light weak transitions. We 
also evaluate a few strong coupling constants appearing in chiral effective lagrangians for heavy and 
light mesons. 


I. INTRODUCTION 

The impressive experimental programme for the study of B decays carried out in recent years has improved our 
knowledge of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and CP violations. In the next few years still more 
abundant data is to come, especially from the dedicated B-factories Belle [jl] and BaBar Q. 

One of the most important goals in these investigations will be a more precise determination of the CKM matrix 
element V u b and, to this end, both exclusive and inclusive b —> u semileptonic transitions will be used. The two methods 
have their own uncertainties. Using the inclusive reaction implies the need to use perturbative QCD methods in the 
region near the end-point of the lepton spectrum, where many resonances are present and perturbative methods 
are less reliable. This difficulty can be avoided by considering exclusive channels and summing them up or taking 
them separately; however the use of the exclusive channels forces to evaluate the hadronic matrix elements by non 
perturbative methods that are either approximate or model-dependent. Examples of these approximations are given 
by non perturbative methods derived from Quantum-Chromo-Dynamics (QCD) first principles, i.e. Lattice QCD and 
QCD Sum rules. The drawback of these methods is the difficulty to improve the precision and to evaluate reliably the 
theoretical error, which follows from the nature of the approximations, i.e. the quenched approximation in Lattice 
QCD and a truncated Operator Product Expansion in QCD Sum rules. Although less fundamental, other approaches 
can be reliably used to give estimates of the hadronic matrix elements that appear in exclusive b —> u transitions and 
we refer here to constituent quark models. At the present stage of our understanding of the hadronic interactions 
from first principles, they offer in our opinion a viable alternative, and the model dependence, which is obvious in 
this approach, can be used to estimate the overall theoretical uncertainty. In this paper we shall follow this route and 
use a Constituent Quark Meson (CQM) model, introduced and defined in || (hereafter referred to as I) to compute 
two semileptonic exclusive heavy to light decays, viz. B to the light vector mesons p and a\. The first decay has 
been recently observed by the CLEO collaboration [|J (see also ||) that has measured the branching ratio of the 
semileptonic decay B —> plv. 


B(B° -> p~tv) = (2.5 ±0.4^ ±0.5) x 10" 4 . (1) 

This decay will be used as a test of the model of Ref. since we do not introduce here any new parameter. On the 
other hand, B —> a\tv is a prediction of this model yet to be tested by the experiment. 

The CQM model developed in I tries to conjugate the simplicity of an effective lagrangian encoding the symmetries 
of the problem together with some dynamical information coming from the observed facts of confinement and chiral 
symmetry breaking. In spite of the simple way in which the dynamics is implemented, the results found in I are 
encouraging. We discussed there the following issues: leptonic constants for heavy mesons, Isgur-Wise form factors 
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for heavy mesons of both negative and positive parity, strong coupling constants of heavy mesons and pions, radiative 
decays of D*, B*\ the comparison with data was satisfactory whenever it was available. 

The plan of the present paper is as follows. The model is briefly reviewed in Section 2. In Section 3 we compute the 
direct contribution to the form factors for the B —> p(is, B —> a\£v decays, i.e. the contribution arising from diagrams 
where the weak current directly interacts with the quarks belonging to heavy and light mesons. In Section 4 we 
compute the strong coupling constants of p and a\ to heavy mesons: these couplings are relevant for the calculation 
of the polar diagrams, i.e. the diagrams where the weak current couples to B and p (or a\) via an intermediate 
particle. These contributions (called polar contributions) are described in Section 5. In Section 6 we present our 
results and compare them with other approaches and with available data. In Section 7 we draw our conclusions and 
in the Appendix we collect formulas and integrals used in the calculations. 


II. THE CQM MODEL 

We begin with a short summary of the CQM model; for a more detailed treatment see I. The model is an effective 
field theory containing a quark-meson lagrangian: 


C = Lit + C-ht ■ 


( 2 ) 


The first term involves only the light degrees of freedom (a constituent quark model for light quarks and mesons was 
originally suggested by Georgi and Manohar ||). To the fields considered in I, i.e. the light quark field x and the 
pseudo-scalar SU (3) octet of mesons 7r, we add the vector meson and axial vector meson octets p p and a p . Considering 
only the kinetic part of the light quarks and mesons as well as the quark-meson interactions at the lowest order, we 
have for Cu- 


C U = + 'tr[F(p) W"] + 

° *9v ^9 a 

+ x{ iDfl hi - ™)x 

+ x(KA* l 'y fl j5 - ihpp^^h, - ih a a> J -'y IJ ,j 5 )x- 


( 3 ) 


Let us discuss the various terms in this equation. The first three terms refer to pions, light vector and axial vector 
respectively. We have £ = exp(*7r// 7r ), E = £ 2 , = 132 MeV; the p and ai field strengths are given by 


A\x)fj.v — dpXjy difXp T , 

where, consistently with the notations of [Q (see also |8)), we write 

m 


n _ • 9v . 
Pv ~ 1 


g v = -f- ~ 5.8 . 

Jit 


By analogy we also write ( m a ~ 1.26 GeV is axial vector meson mass): 


■ 9A . 

V2 


m a 


9A = -f- ^ 9.5 . 
J 7T 


( 4 ) 


( 5 ) 


( 6 ) 


Here p, a are hermitean 3x3 matrices of the negative and positive parity light vector mesons. The fourth term in 
Eq. (||) contains the light quarks, with D M = <9 M — iVp and 




( 7 ) 


Therefore it gives both the kinetic term of the light quark and its coupling to an even number of pions. For m , in I 
we took the value m = 300 MeV (for non-strange quarks). The last three terms describe further interactions between 
light quarks and light mesons. The coupling of the light quark to an odd number of pions is mediated by 




( 8 ) 


Moreover, consistently with a low energy theorem for pions, we put = 1. Concerning the interactions of vector 
particles, we put h„ = 


\[2 

9vf P 


i ha — 


V^rric 

9Afa 


where f p and f a are the leptonic constants. For the p leptonic constant we 
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use f p = 0.152 GeV 2 , as given by p°, u> decay into e + e~. For f a a phenomenological determination using r —> v T nnn 
was obtained in |J, i.e. f a = 0.25 ± 0.02 GeV 2 , a result which agrees with the one found by QCD sum rules |n|. On 
the other hand from lattice QCD one obtains f a = 0.30 ± 0.03GeV 2 jn|]. Since 1 / f a multiplies all the amplitudes 
involving the a± meson, the uncertainty on f a will induce a normalization uncertainty on all the amplitudes involving 
the light axial-vector meson. We note that our choice for h p and h a implements the hypothesis of the Vector and 
Axial-Vector Meson Dominance. Numerically we find: 

h p ~h a ~ 0.95 . (9) 

We also observe that our choice for the normalization of the light axial vector meson field Eq.(|^) is conventional since 
gA disappears from the physical quantities (in Q gA = gv is assumed). We also differ from the phenomenological 
analyses of Ref. ||] since we do not assume the current algebra relations to^ = 2to 2 and f a = f p that seem to have 
substantial violations. 

Let us now discuss Che, he. the part of the lagrangian that contains both light and heavy degrees of freedom, in 
particular the heavy quark ( Q ) and mesons ( Qq ). According to the Heavy Quark Effective Theory (HQET) |Q, in 
the limit toq —> oo, these mesons can be organized in spin-parity multiplets. We shall consider here the negative 
parity spin doublet (P, P*) (e.g. B and B*) and its extension to P-waves, i.e. the doublet containing the 0 + and 1 + 
degenerate states (Po,P*'). Incidentally, we note that HQET predicts another low-lying multiplet, comprising two 
degenerate states with 1 + and 2 + |l3|], which is of no interest here. In matrix notations these fields can be represented 
by two 4x4 Dirac matrices H and S , with one spinor index for the heavy quark and the other for the light degrees 
of freedom. An explicit matrix representation is, for negative parity states: 

h = [p;r - P75] (io) 

(.H = 7 oP^ 7 o), whereas, for positive parity states: 

s = ^[pr;7 M 75-p 0 ] ■ (u) 

In these equations v is the heavy meson velocity, tPP* = t !, 'P( ; ( = 0; P*^, P, P£' and Pq are annihilation operators 
normalized as follows: 

<0|P|Qg(0-)> = \/Mh (12) 

(0\P^\Qq(l-)) = e»y/M^ , (13) 

with similar equations for the positive parity states (Mp = Mp = Mp* is the common mass in the H multiplet). 
With these notations the heavy-light interaction lagrangian is written as follows: 

Chi = Qviv ■ dQ v - (x(H + S)Q V + h.c.) 

+ + S)(H - S)) , (14) 

where Q v is the effective heavy quark field of HQET and we have assumed that the fields H and S have the same 
coupling to the quarks, which is a dynamical assumption based on a simplicity criterion (in I we justify it on the basis 
of a four quark Nambu-Jona Lasinio interaction by partial bosonization ]lj[|). After renormalization of the heavy 
fields H and S 0 one obtains the kinetic part of the heavy meson lagrangian in a form that is standard for heavy 
meson effective chiral theories 0 : 

Che = Tr H{iv ■ d — A p)H + Tr S(iv ■ d — A s)S . (15) 

Here Ah and As are the mass difference between the meson and the heavy quark at the lowest order; typical values 
considered in I are Ap = 0.3 — 0.5 GeV. Ap and A 5 are related: for example, for A# = 0.4 GeV one obtains the 
value As = 0.590 GeV ||. These values correspond to a value for the 5—multiplet mass m = 2165 ± 50 MeV; these 
states, called in the literature (for the charmed case) Dq, D* r , have not been observed yet, since they are expected to 
be rather broad. H and S are the renormalized fields and are given in terms of the bare fields H, S by 


H= * 

V Zp 

(16) 

s= 

V^s 

(17) 
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Zh, Zs are renormalization constants that have been computed in || with the results (the integral I 3 can be found 
in the Appendix): 


Zh 


Z 


-1 

s 


(A H + m) dh ^ H) +J 3 ( A h ) 
ol\ h 

{A s -m)^^ + h(A s ) 


(18) 

(19) 


where m is the constituent light quark mass. 

Let us finally discuss the way to compute the quark-meson loops arising from the previous lagrangian. As we have 
seen, the CQM model describes the interactions in terms of effective vertices between a light quark, a heavy quark 
and a heavy meson (Eq. <ED- We describe the heavy quarks and heavy mesons consistently with HQET; for example 
the heavy quark propagator is given by 


i 

v ■ k + A 


( 20 ) 


where A is the difference between the heavy meson and heavy quark mass and k is the residual momentum arising 
from the interaction with the light degrees of freedom. 

The light quark momentum is equal to the integrated loop momentum. It is therefore natural to assume an 
ultraviolet cut-off on the loop momentum of the order of the scale at which the chiral symmetry is broken, i.e. A ~ 1 
GeV (in I we assumed the value A = 1.25 GeV). Since the residual momentum of the heavy quark does not exceed 
few units of A qcd in the effective theory, imposing such a cut-off does not cut any significant range of ultraviolet 
frequencies. We also observe that the value of the ultraviolet cut-off A is independent of the heavy quark mass, since 
it does not appear in the effective lagrangian. 

Concerning the infrared behavior, the model is not confining and thus its range of validity cannot be extended below 
energies of the order of A qcd- In order to drop the unknown confinement part of the quark interaction one introduces 
an infrared cut-off p. These parameters appear in the regularized amplitudes; as discussed in || (see also Q) we 
have chosen a proper time regularization; the regularized form for the light quark propagator (including integration 
over momenta) is 


[ d 4 k E , 2 1 , -> [ d 4 k E / V#1 dse- s ^ +m ^ , 

J k E + m 2 J J 1/A 2 


( 21 ) 


where /.1 and A are the infrared and ultraviolet cut-offs. For /i in I we assumed the value: p = 300 MeV. For a different 
choice of the cut-off prescription in related models see [^5| , jl(J . 


III. B -> p AND B -*• Ai FORM FACTORS: EVALUATION OF THE DIRECT CONTRIBUTIONS. 

The form factors for the semileptonic decays B —> p£u and B —> a\iv can be written as follows (q = p — p'): 

< P + {tWiP') |m7/x(1 - ls)b\B°(p) > = -^y^-e F a/J£’V/ 


where 


ms + m p 
- ie* (m B + m p )A 1 (q 2 ) 


+ «(e* • q) 


{r + p,) - a 2 (,’) 


ms + m p 


2 m 


+ *(e* ’ ?)—/9/4^3Or) - A o (9 )] 


3 / 2\ m B + m p a / 2\ m B m p 

Mq ) = — Mq ) — Mq ) 


2 m 


2 m n 


( 22 ) 


(23) 


and 
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< af(e(A),p , )l9 , 7^( 1 - l5)b\B(p) > = 

m B + m a 

- *e*(m B + m a )V 1 (q 2 ) 

+ i(e*-q) (P+ / )/J F 2 ( g 2 ) (24) 

m B + m a 

+ i{e* ■ q)^^q tl [V 3 (q 2 ) - V 0 (q 2 )] , 

where m a is the a± mass and 


t / , m B + m a T/ , 2 , m B - m a T ^ , 2 , 

niq > = ~^r Vi{q > - ~^r VAq > • 

We note that, having used this parameterization for the weak matrix elements |l7|| , at q 2 
hold 


(25) 


0 the following conditions 


MO) = 4>(0) ( 26 ) 

v^(0) = Vo(0) ■ ( 27 ) 


The contribution we consider in this section arises from diagrams where the weak current couples directly to the 
quarks belonging to the light and heavy mesons (see Fig. 1). 

These diagrams are computed using the rules described in the previous Section. The results of a straightforward, 
but lengthy calculation are as follows: 


where 


and 


f P 

2 m 2 


Zh 


V (q ) = — -J— \ -(fli - mZ) ( m B + m p ) 


m B 


1 


A i(q 2 ) = —r^-VZHms- 

}p m B + m p 

[(m 2 + mm p u))Z — umi p Q \ — m p fl 2 — 2n 3 — 

II4 — O5 — 2tjflg] 

fig 


Mia 2 ) = ^ 

jp 


Z B 

m B 


mZ — fli — 2- 


{m B + m p ) 


A°(q 2 ) = - y-\fZ H m B 

i P 


q r i \ 

fli ( m p Co + 2 m — 5 -+ mpVL 2 + 

m B m B J 

r\ 


2f2 3 T ST 4 (1 — 2—2~) T fig T 2flg I cu — 


ry, 2 r l 

Z ( m — m - 

m B 


q B 

mm p Cu) 


m B m p 


wU w p- g 2 
2 m B m p 

2 2 2 
m 2 B — q 2 — m 2 

2 


(28) 


(29) 

(30) 


(31) 


(32) 


(33) 


and the functions Z , fl, are given by the formulae of the Appendix with Ai = A B , A 2 = Ai — m p u>, x = to p ; m is 
the constituent light quark mass. 
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The calculation for the B —> a± transition is similar. The results are: 

A D {q 2 ) = (^i ~ mZ - —12 2 ) ( m B + m a ) 

Ja V rn B \ rn a ) 

1 


where now: 


y D{q 2 ) = --- 

Ja rn B + Wo 

[(—m 2 + mm a Lu)Z + 2mfli — um a fl 1 + 
+(2mu — m a )fl2 — 2 f 23 — 114 — ^5 — 2coV.e] 

V°(q 2 ) = (mZ — Qi — 2 — + 2—^2) {m B 

Ja V tub \ m a m a J 


+ rn a ) 


V 0 D (q 2 ) = — j-y/Z H m B 

Ja 


fii ( m a uj + 2 m-J, --- 2 m 


„/ 


1 ± 

m B 


1^2 (TIT a -f- 2 TYl. - — 2?TIUj) T 2123 T 12 4 (\ — 2——~\~ 11") — 

V m B m a J v m B J 


212 , 


6 W 


m B m a 


u> = 


+ Z ( m 2 + m—i mm a ui 

m B 


m 2 B + m 2 — q 2 

2 m B m a 
2 2 2 
m s - r - m z a 


( 34 ) 


(35) 

(36) 


(37) 


(38) 

(39) 


The previous results for the form factors can be used directly for the numerical analysis. In order to allow an easier 
way of using our results we have fitted these formulas by the simple parameterization: 


F D (q 2 ) 


F D { 0) 



(40) 


for a generic form factor F D (q 2 ); of ,&f have been fitted by a numerical analysis performed up to q 2 = 16 GeV 2 , 
both for p and ai mesons. We have collected the fitted values in Table [j] We note explicitly that the results for 
B —> a± form factors at q 2 = 0 are proportional to the factor (0.25 GeV 2 // a ). Besides the normalization uncertainty 
due to f a , we estimate a theoretical error of 15% on these parameters. 


IV. STRONG COUPLING CONSTANTS 

In this Section we compute the strong couplings FJFJp, HSp , HFJai, FJSa\. As discussed in the introduction they 
are relevant for the calculation of the polar contribution to the form factors. We parameterize these couplings by 
considering the following effective Lagrangians (we follow the notations introduced in Q|): 

Chh p = i\Tr(HHa^F(p)^) - ipT\{HH^ Plx ) (41) 

Chs p = -<Tr (SH^Pp) + ipTr(SHa^F(p)^) (42) 

C HHai = -iC A Tr(HH'fa li ) + i p A Tr (H Ha^F (a) (43) 

C HSai = i\ A Tr(SHa^F(a)^) - 0 A Tr^H'fa ft ) . (44) 

The strong couplings can be computed by matching the effective meson lagrangian of Eqs.(^)-(|4j) with the quark- 
meson lagrangian ©> i-e- considering triangular quark loops with external legs representing light and heavy mesons. 
The calculation is similar to the one of the previous section. The results are as follows: 

nr 2 

A = p Z H (~ni + mZ) (45) 

v 2 gvfp 

2 

Tn 

/3 = V2 — ?-Z H [ 2 m£h + m p n 2 + 212 3 -12 4 +12 5 - m 2 Z] . (46) 

9v Jp 
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Here the functions Z, Qj are given by the formulae of the Appendix with Ai = A h, x = m p , lo = m p /{2ms) (we 
keep here the first 1 /toq correction). Moreover 


H = l np , ^Z H Z s {-n i - 2^. + mZ) 

V^gvfp m p 

fX 2 

C = - ~r~ \JZ\i Zg (m p O ,2 ± 2 H 3 ± H 4 ± H 5 — m~Z) . 

gvfp 

Here the functions Z, flj are given by the formulae of the Appendix with Ai = A h, A 2 = A 5 , x 
u> = (Ai — A 2 ) jTrip. For the axial-vector ai couplings to H and S states we find 

t y/Z H Z s {-Sh + 2H 2 — + mZ) 

V2g A fa rn a 

2 

V2 — y- \JZHZs{m a £l2 ± 2 H 3 — H 4 + H 5 ± to 2 Z) , 
gAja 

where Z , Qj are given by the formulae of the Appendix with Ai = Ah, A 2 = A 5 , x = m a and w = (Ai 
Moreover 


Aa = 
Pa = 


(47) 

(48) 
= m p and 


(49) 

(50) 
- A 2 )/m a . 


HA 

U 


V2g A f< 

V2m 2 a 


Z H to [ Z + 2 


no 


-n ± -2- 




9A fa 


Zh { —2777/^1 + 777-0,^2 2^3 ^4 ^5 777 Z ) , 


where Ai = A#, x = m a: cu = m a / (2ms) Numerically we get the following results: 


A = 

0.60 GeV -1 

Aa 

P = 

-0.86 

Pa 

g = 

0.16 GeV" 1 

HA 

C = 

0.01 

Ca 


0.85 x (0.25 GeV 2 /f a ) GeV -1 
-0.81 x (0.25 GeV 2 /f a ) 

0.23 x (0.25 GeV 2 /f a ) GeV” 1 
0.15 x (0.25 GeV 2 /f a ) . 


(51) 

(52) 


A discussion about the theoretical uncertainties of these results is in order. We have explicitly written down the 
dependence on f a of the strong coupling constants involving the light axial-vector meson, since as noted before, this 
is a major source of theoretical uncertainty for these constants. Another source of spreading in the reported values 
is the variation of A h in the range A h = 0.3 — 0.5 GeV (we use A h = 0.4 GeV in the calculation). This produces 
a significant uncertainty only for C, Pa-, Cx since we obtain f = 0.01 ± 0.19, Pa = — 0 . 81 ^q ; 2 4 and £ a = 0.15lg ;[| 
while in the other cases only a few percent variation is observed. For the other constants: A, g, Aa, fJ-A, a theoretical 
uncertainty of ±15% can be guessed. This estimate follows for example from a different evaluation of the A parameter 
performed in I. For other determinations of the coupling constant A see jn| (QCD sum rules and light cone sum rules) 
and ||l9|l (this paper uses data from D* decays together with vector meson dominance). 


V. B -> p AND B ->• Ai FORM FACTORS: EVALUATION OF THE POLAR CONTRIBUTIONS. 

The polar contributions are given by diagrams where the weak current is coupled to B and to the light vector or 
axial vector meson by an intermediate heavy meson state (see Fig. 2). These diagrams, because of the heavy meson 
propagator, produce a typical polar behavior of the form factors, in the form 

= (53) 

m'p 

This behaviour is certainly valid near the pole; we assume its validity for the whole q 2 range, which can be justified 
on the basis of the minor numerical role of this polar contribution, as compared to the direct one, in the region of low 
q 2 , at least for the form factors Af,A^ (see the numerical results at the end of this Section). The assumption ( ]TS| ) 
cannot be made for the form factors A //( q 2 ) and V 0 P ( q 2 ), as we discuss below, and is also less reliable for A p (q 2 ) and 
V p (q 2 ) (see Table 0). 
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The values at q 2 = 0 in Eq. (^) can be easily computed in term of the strong coupling constants defined in the 
previous Section and using the leptonic decay constants F and F + that give the coupling of the intermediate states 
to the currents. Neglecting logarithmic corrections, F and F + are related, in the infinite heavy quark mass limit, to 
the leptonic decay constant fs and f + defined by 


{^\qi^lbb\B(p)) = ip p f B (54) 

(0\rfb\B o (p)) = p p f + , (55) 


by the relations fs = F / ^/m B and f + = F + / ^/TOb 0 (B 0 is the S state with J p = 0 + and bq content). These 
couplings have been computed in |j| with the results given in Table Ill for different values of the parameters. 

For the values F p ( 0) we obtain the following results for the B —> p transition. 


F p (0) 

Af(0) 

Af( 0 ) 


— V2gv^F 


mB + m p 
3/2 


y/2 m B gvF + . 

- - ---r(C - 2 lIUJUlp) 

m B o ( m B + m p ) 

m B 0 


(56) 

(57) 

(58) 


where u> = ms /(2m p ). For A p ( q 2 ), we have to implement the condition contained in Eq. (^6[); for instance a possible 
choice is 


A p (q 2 ) = A p (0)+g v (3F 


1 


m p y/2m B m 2 B - q 2 


For the B —> a\ transition we have: 


^ P ( 0 ) 

Vf(0) 


3/2 


-V 2 g A X A F +mB+m ' 

y/2m B 9AF 


mB(mB + m a ) 

-Aig^P^sFUL 


(C A - 2fJ, A UJma) 


(59) 


(60) 

(61) 

(62) 


where u> = mB /(2m a ). Similarly to the previous discussion for A p (q 2 ) we can put for instance: 


V p (q 2 ) = V p (0)+g A p A F + 


m a y/2m B m 2 Bo - q 2 


(63) 


We note that Eqs. (59) and (|63|) have been written down only as an example of a possible behaviour of these form 
factors satisfying the given constraints. For massless leptons they do not contribute to the semileptonic width and 
can be neglected. 

Numerically we obtain the results in Table [H] where we have also reported the values of the pole masses for all the 
form factors except A p (q 2 ) and V p (q 2 ) because of Eqs. ( |59| ) and (^3|). Similarly to the previous analyses an overall 
uncertainty of ±15% can be assumed. In Fig. 3 we plot the form factors A\ and A 2 for the semileptonic decay B —» p. 
In Fig. 4 are shown the form factors A, Vi and V 2 for the semileptonic decay B —> a±. Since the behaviour in Eqs. (|5S| ) 
and ( p3|) is only guessed, we have not included the form factors A p ( q 2 ) and V p ( q 2 ) in Fig. 3 and 4; in addition V{q 2 ) 
is not reported in Fig. 3 since our prediction is affected by a large error (see the discussion in the next section). Note 
that the theoretical error is not included in Fig. 3 and 4; one should refer to the numbers in Tables | and ||. 


VI. BRANCHING RATIOS AND WIDTHS 

In this Section we compute the branching ratios and widths for semileptonic decays using the numerical results for 
form factors reported in Table I,II. Let us first compare our results for the B —> p form factors with those obtained 
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by other methods (see Table IV). These form factors (as well as those concerning the transition B —> ai) are obtained 
by adding the direct and polar contributions: 

F(q 2 )=F D (q 2 )+F p (q 2 ), (64) 

where F D (q 2 ) was introduced in Section III and F p (q 2 ) in Section V. Our result for the vector form factor V p (q 2 ) is 
affected by a large error since it arises from the sum of two terms opposite in sign and almost equal in absolute value. 
Apart from this large uncertainty, our results are in relative good agreement with the results of QCD sum rules, but 
they are in general higher than those obtained by other approaches. For the B —> plv decay width and branching 
ratio we obtain (using V u b = 0.0032, tb = 1-56 x 10 -12 s): 

B(B° -> p+tv) = (2.5 ± 0.8) x 10~ 4 

T 0 (B° -► p + £v) = (4.4 ± 1.3) x 10 7 s" 1 

T+(I3 0 -> p+£v) = (7.1 ± 4.5) x 10 7 s" 1 

r_(5° -► p+iu) = (5.5 ± 3.7) x 10 7 s" 1 

(T+ ± r_)(5° -> p+tv) = (1.26 ± 0.38) x 10 8 s" 1 (65) 

where To, r+, T_ refer to the p helicities. The branching ratio for B —► p£i/ is in agreement with the experimental 
result quoted in the introduction, Eq. ([]]). 

Let us now discuss the theoretical uncertainty of these results. The large error of V p (0) affects significantly the 
values of T + and T_, whose errors are correlated; it has however no effect on To and a small effect on the branching 
ratio, which increases at most by 8%. The theoretical uncertainties on A p (0) and A^(0) are likely to be related. To 
get the theoretical error on the widths we have added in quadrature the error induced by V p (0) and a common ±15% 
error on A p (0) and A%( 0). 

Having used the decay B —> piv as a test of the CQM model, we can now consider the B —> a\tv semileptonic 
decay. The results obtained are: 


GS 

toi 

o 

1 

-> a+lv) = (8.4 ± 1.6) x 10 “ 4 


r o(B° - 

-4 a+tv) = (4.0 ±0.7) x 10 s s " 1 


r +(B° - 

a\£v) = (4.6 ±0.9) x 10 7 s " 1 


r_(H° - 

a+£v) = (0.98 ± 0.18) x 10 8 s " 1 

( 66 ) 


where To, T + , T_ refer to the a\ helicities. We have included in the determination of these decay widths only the 
normalization uncertainty arising from f a ; the lower values correspond to f a = 0.30 GeV 2 while the higher values to 
f a = 0.25 GeV 2 . One should also take into account the theoretical errors arising from the values of the form factors 
at q 2 = 0; they are more difficult to estimate reliably and are not included here. In any case the over-all theoretical 
uncertainty is larger (presumably by a factor of two) than the one reported in the previous formula. 


VII. CONCLUSIONS 

The main conclusion of this paper can be read from eqs.(| 66 |). We predict a branching ratio for the decay B° —> a p £v 
significantly larger than the branching ratio for B° —> p + £v\ in spite of the theoretical uncertainties inherent to the 
CQM model, which we have discussed in the previous Sections, this is a remarkable outcome. A consequence of this 
result is that the B —> a\£v decay channel might account for around 50% of the semileptonic B —> X u £v decay channel 
(evaluated, for example, by the parton model), whereas the B —> plv decay channel adds another 15%; given the 
relevance of these results for the determination of V u b, it would be interesting to test these predictions in the future 
by other theoretical methods and, hopefully, by experimental data. 


VIII. APPENDIX 

In the paper we have introduced several integrals and parameters that we list in this Appendix. 
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d 4 k 


N r 


(v ■ k + A + ie) 
^ ds 


L % (i + ™ 2 - a2 ) [1+erf(a ^ )1 


16 7r 3 / 2 J!/ A 2 s 3 / 2 
A N c m 2 m 2 to 2 

^ - " " 4 V A 2 >. .2 / 


107T 2 


A 2 ’ n 2 


h = 
I'i = 


iN c l' reg d 4 k N c m 2 . m 2 m 2 
(k 2 - m 2 ) = 16^- (_1) A^V 


167T 4 

iA r c 


re g 7.2 

d 4 k 


2 2 
m m 


h = - 

H A) = - 


167T 4 J 

iN r r s d 4 k 


_= ^! r (-2 — —) 

(fc 2 - m 2 ) 8 tt 2 v ’ A 2 ’ /i 2 


N r 


167T 4 

iN r . 


, 7/6 7/6 . 

(fc 2 -m 2 ) 2 ~ 16 tt 2 ' ( °’ A 2 ’’^ 2 " ) 
f reg 


107T 4 7 (fc 2 — m?){v ■ k + A + ie) 


N r 


r i/m 2 


16 7T 3 / 2 Ji/A 2 S 3 / 2 


ds e _ s(m 2_ A 2 } 


(l + erf(Av / s)) 


j/V f reg 


d 4 fc 


Nr 


( k 2 — m 2 ) 2 {v ■ k + A + ie) 
ds 


[1 + erf (Ay's)] . 


167T 3 / 2 7i/ A 2 s 1 / 2 

where T is the generalized incomplete gamma function and erf is the error function. Moreover, having defined: 

Ai (1 — x) + A 2 x 


we have: 


cr(x, A 1 i A 2 ,oj) = 


I 5 ( Al ,A 2 , u ) = ^f 


= dx 


i/l + 2 (u> — 1) x + 2 (1 — u>) x 2 

g _ d^k _ 

(fc 2 — m 2 )(v ■ k + Ai + ie)(v' ■ k + A 2 + ie) 
1 


(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 

(73) 


1 + 2a; 2 (l — u) + 2x(u — 1) 


6 


r i/m 2 


L167T 3 / 2 y i/A 2 

6 r 1 /^ 


167T 2 


ds a e s ( m a ^ s l ^ 2 (1 + erf(crv / s)) + 
ds e ~ s (™ 2 - 2 ° 2 ) s - 4 . 


'1/A 2 


(74) 


We also define, if q M = xv'^, u = v ■ v', A 2 = A 1 -iw, the formula: 

„ iN c r s d 4 k 

z = 


167T 4 J ( k 2 — m 2 )[{k + q) 2 — m 2 ](v ■ k + Ai + ie) 

_ h(Ai,x/2,u>) - I 5 {A 2 ,-x/2,u>) 

2x 

We use in the text the following combinations of the previous integrals: 

K 1 = m 2 Z-I 3 {A 2 ) 

K 2 = A \Z - ^ JV "'“’ [io x + 2Ai] 


( 75 ) 


K 3 = — Z + 


k 2 ry h( x/2) ~ I 3 (~ x/2) f 
4x 

X“ „ . / 3 (Ai) - 3 / 3 (A 2 ) uj 


+ [^1 d 3 (A x ) — A 2 / 3 (A 2 )] 
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K 4 = 

Q 1 — 
f^2 = 

f^3 = 


n 4 


1 1 — 


fig = 


sAt A 1 [J 3 (A 1 )-/ 3 (A 2 )] h[x/2)-h{-x/2) 

2 2x 4 

h(—%/2) - I 3 (x/2) +u;[/ 3 (A 1 ) - J 3 (A 2 )] [A 3 - a;a;/2]Z 

2x(l — w 2 ) 1-w 2 

-/ 3 (Ai) + / 3 (A 2 ) - w[/ 3 (-s/2) - h{x/2)} [x/2 - A 4 u\Z 


2x(l — u 2 ) 

Ki 2 ujK a -K 2 - K :i 


1 — UJ 2 


2 2(1 -uj 2 ) 

-K a , 3A' 2 - GujK 4 + K 3 {2ui 2 + 1) 
2(1 -w 2 ) + 2(1 — w 2 ) 2 

—-Ki , 3if 3 -6wif 4 + ^2(2w 2 + l) 
2(1 — u > 2 ) + 2(1 — w 2 ) 2 

K 4 uj | 2/\ 4 (2w 2 + 1) - 3w(AT 2 + AT 3 ) 
2(1 — w 2 ) + 2(1 — to 2 ) 2 


(76) 
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Fig. 3 - Form factors A\ (continuous line) and A 2 (dashed line) for the semileptonic decay B — > plv. 



Fig. 4 - Form factors A (continuous line), V\ (dashed line) and V 2 (dot-dashed line) for the semileptonic decay 
B —> a\lv. 
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V D 




A u 

vf 

V? 

V u 

v o 

F u (0) 

0.83 

0.69 

0.81 

0.33 

1.62 

1.13 

1.13 

1.13 

dp 

0.93 

0 

0.87 

2.9 

1.13 

0.18 

1.3 

1.9 

b F 

0.02 

0 

—0.17 

2.6 

0.12 

0.04 

3.8 

0.93 


TABLE I. Parameters of the direct contribution to the various B form factors for B —> p and B —> ai decays. The values 
F d ( 0) for the B —» ai transition (last four columns) should be multiplied by the factor 0.25 G eV 2 // a . The theoretical 
uncertainty is ±15%. 



v p 

< 

A? 


A p 

vf 

vf 

V p 

F^(0) 

mip 

-0.84 

5.3 

-0.11 

5.5 

-0.15 

5.5 

-0.019 

-1.48 

5.5 

-0.32 

5.3 

-0.57 

5.3 

0.07 


TABLE II. Parameters of the polar contribution to the various B form factors for B —> p and B —> ai decays. Pole masses 
in GeV. The values F p ( 0) for the B —> ai transition (last four columns) should be multiplied by the factor (0.25 GeV 2 // a ). 
The theoretical uncertainty is ±15%. 


A h 

F 

F+ 

0.3 

0.33 

0.22 

0.4 

0.34 

0.24 

0.5 

0.37 

0.27 


TABLE III. F and F + for various values of Ah- Ah in GeV, leptonic constants in GeV 3 / 2 . 



This work 

Potential Model [^0| 

LCSR @ 

SR@ 

Latt. + LCSR @ 

V p {Q) 
AK 0) 
^(0) 
4J(0) 

-0.01 ±0.25 

0.58 ±0.10 

0.66 ±0.12 

0.33 ±0.05 

0.45 ± 0.11 

0.27 ± 0.06 

0.26 ± 0.05 

0.29 ± 0.09 

0.34 ± 0.05 

0.26 ± 0.04 

0.22 ± 0.03 

0.6 ± 0.2 

0.5 ± 0.1 

0.4 ± 0.2 

0.24 ±0.02 

u.oo_ 0 .05 
r» 97+0.05 
u - z ‘-0.04 

0 26 +0 ' 05 

U.ZO_ 0 .03 

0 SO+ 0 - 06 

U.OU_ 0 .o4 


TABLE IV. Form factors for the transition B —> p at q 2 = 0. Our results are compared with the outcome of other theoretical 
calculations: potential models, light cone sum rules (LCSR), QCD sum rules (SR), calculations involving both lattice and 
light cone sum rules. The large error of V p (0) in our approach is due to the large cancellation between the direct and polar 
contribution. 
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